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We present an update of our analysis of short-baseline neutrino oscillation data in the framework 
of 3+1 neutrino mixing taking into account the recent update of MiniBooNE antineutrino data and 
the recent results of the MINOS search for i/^ disappearance into sterile neutrinos (the more com- 
plicated 3-1-2 neutrino mixing is not needed since the CP-violating difference between MiniBooNE 
neutrino and antineutrino data has diminished) . The results of our fits of short-baseline neutrino os- 
cillation data including the MiniBooNE low-energy anomaly (now present both in the neutrino and 
antineutrino data) leads to a strong tension between appearance and disappearance data. Hence, 

it seems likely that the low-energy anomaly is not due toV^l — J-Ve transitions. Excluding the Mini- 
BooNE low-energy anomaly, appearance and disappearance data are marginally compatible. The 
global analysis has the best-fit point at Am^i « 5.6 eV^, which is rather large in comparison with 
cosmological bounds, but there are three regions within la at Amli ^ 1.6, 1.2, 0.91 eV^. We also 
show that the data on the Gallium neutrino anomaly favor values of Am|i larger than about 1 eV^. 

PACS numbers: 14.60.Pq, 14.60.Lm, 14.60.St 



I. INTRODUCTION 

The possible existence of sterile neutrinos is an ex- 
citing possibility which could open a powerful window 
on our view of the physics beyond the Standard Model. 
The short-baseline (SBL) neutrino oscillation experiment 
LSND [1] discovered in the late 90's a signal which 
can be due to — oscillations generated by a 
neutrino squared-mass splitting Arngg^ of the order of 
1 eV^, which is much larger than the well-established so- 
lar (SOL) and atmospheric (ATM) squared-mass split- 
tings, Am|oL = (7.6±0.2) x lO^^eV^ [2] and Am^^M = 
2.32^0,11 X 10"^ eV^ Q- In order to have more than two 
independent squared-mass splittings the number of mas- 
sive neutrinos must be larger than three. In this case 
the flavor neutrino basis is composed by the three known 
active neutrinos Ve, f/i, i^t and by one or more sterile 
neutrinos i/si, ^^2, ■ ■ ■ , which do not have standard weak 
interactions and do not contribute to the number of ac- 
tive neutrinos determined by LEP experiments through 
the measurement of the invisible width of the Z boson, 
Na = 2.9840 ± 0.0082 

Schemes of neutrino mixing with sterile neutrinos have 
been studied by several authors (see Refs. 0-111), with 
more attention to the simple schemes with one or two 
sterile neutrinos (four- and five-neutrino mixing, respec- 
tively). Since the three active neutrinos must have large 
mixing with the three massive neutrinos which generate 
AmgQL and Am\rj,yi and no effect of sterile neutrinos 
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has been seen in solar and atmospheric neutrino data, 
the mixing schemes with sterile neutrinos must be per- 
turbations of the standard three-neutrino mixing scheme 
in which the three active neutrinos Ve, v^i, Vt are superpo- 
sitions of three massive neutrinos vi, vi, vz with respec- 
tive masses toi, m^, m^, such that ArrigQ^ 



'ATM — l^"^31l — I^"''32|J "i"^^ "'"fej ~ '"fc ■■"J- 

Moreover, standard analyses of the Cosmic Microwave 
Background and Large-Scale Structures data constrain 
the neutrino masses in the case of three-neutrino mixing 
to be much smaller than 1 eV and are compatible 

with the existence of one or two sterile neutrinos which 
have been thermalized in the early Universe |T3| only 
if the masses of the additional, mainly sterile, massive 
neutrinos are smaller than about 1 cV • Also Big- 

Bang Nucleosynthesis data are compatible with the ex- 
istence of sterile neutrinos which have been thermalized 
in the early Universe [l^, [2l| , with the indication how- 
ever that schemes with more than one sterile neutrino 
are disfavored [l^ [2^. Hence, the schemes with ster- 
ile neutrinos which are currently under consideration are 
the 3-1-1 and 3-1-2 schemes in which Ve, i^^i, Vr are mainly 
mixed with v\, V2, I'a, whose masses are much smaller 
than 1 eV and there are one or two additional massive 
neutrinos, and v^, which are mainly sterile and have 
masses of the order of 1 eV. Short-baseline oscillations 
corresponding to the LSND — > f/g signal are generated 
by the large squared-mass differences Ato|]^ and Amgj^. 

The LSND signal was not seen in the KARMEN 
li, NOMAD i^f, ^ i^e ^ and MiniBooNE 
25l | short-baseline experiments. However, in 
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Am2i and 



July 2010 the interest in the LSND signal has been re- 
vived by the observation of a compatible signal in the 
MiniBooNE — > short-baseline experiment f26j . The 
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MiniBooNE and LSND antineutrino data have been an- 
alyzed in several papers in conjunction with the data of 
other short-baseline experiments in the framework of 3-1-1 
and 3-1-2 schemes [27F30]. In this paper we update the 
3-1-1 analysis presented in Ref. [30,] by taking into account 
the update of MiniBooNE antineutrino data presented in 
Refs. HH, IHl and the recent results of the MINOS search 
for Vi, 



disappearance into sterile neutrinos 
In 3+1 schemes we have the squared-mass hierarchy 



(1) 



and we consider four-neutrino mixing as a perturbation 
of three-neutrino mixing: 
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The effective flavor transition and survival probabilities 
in short-baseline experiments are given by 



SBL 



= sin 2i?Q^ sin 



AE 



SBL 



for a, /3 = e, /x, r, s, with 



sin2 2t?,^ =4|;7„4n;7^4|', 
sin^ 2i?„„=4|;7„4|'(l-|C^«4|') 
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In this paper we do not consider 3-1-2 neutrino mixing, 
because the new MiniBooNE antineutrino data 31, 33 
do not show a sufficient difference from the MiniBooNE 
neutrino data [1^ to motivate the consideration of the 
much more complicated 3-1-2 neutrino mixing which 
could explain a difference throu gh C P violation in short- 
baseline oscillations [11 111, [salsJls^l ■ In Ref. [s^l we 
have shown that such an effect would be the main moti- 
vation for preferring 3-1-2 mixing over 3-1-1 mixing. 

The plan of the paper is as follows. In Sections HIl 
and mil we describe, respectively, our analysis of the new 
MiniBooNE antineutrino data and that of MINOS data 
on the search for v^^ disappearance into sterile neutrinos. 
In Section IIVI we present the results of the global fit of 
short-baseline oscillation data and in Section |V] we draw 
our conclusions. 



II. MINIBOONE 

The MiniBooNE collaboration presented recently a 
preliminary update of their antineutrino data obtained 
with 8.58 X 10^0 Protons on Target (PoT) which 
increases the amount of data by a factor of about 3/2 
with respect to the data published in Ref. [1^, which 
were obtained with 5.66 x 10'^° PoT. The new data show 
two interesting new features: 

1. The antineutrino data have an anomalous low- 
energy excess similar to that of the neutrino data 
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FIG. 1. Allowed regions in the sin^ 2'!9efi-Am4i plane ob- 
tained from the fit of MiniBooNE antineutrino data [3ll . [3^ , 
including the low-energy bins from 200 MeV to 475 MeV. The 
best-fit point at sin^ 2-0^^1 = 0.005 and Am|i = 4.68 eV^ is in- 
dicated by a cross. 



2. The — ?> signal in the three energy bins from 
475 MeV to 800 MeV has sHghtly diminished with 
respect to that published in Ref. [H]. 

The first new feature raises the interesting question 
if the MiniBooNE low-energy anomaly can be fitted by 
oscillations. In order to answer to this question we con- 
sider the fit of MiniBooNE data with and without the 
three low-energy bins from 200 MeV to 475 MeV. The 
second new feature may be the consequence of the fluc- 
tuations of the signal around the true value. A conse- 
quence of the new data is that the difference between 
MiniBooNE neutrino and antineutrino data has dimin- 
ished, lessening the need of CP violation that was sug- 
gested by the previous data. Hence, there are less motiva- 
tions for considering 3-1-2 neutrino mixing pTi, [29l. Isol, Is^- 
ISTj . in which v^^ v^. and oscillations can be 

different if \Yi\\Ul^U ehU^^^ ^ 0. Moreover, as we ar- 
gued in Ref. [30,] , 3-1-1 neutrino mixing is preferable on 
3-1-2 mixing for its simplicity and for the natural corre- 
spondence of one new entity (a sterile neutrino) with a 
new effect (short-baseline oscillations). In Ref. jiCi] we 
have also shown that the improvement of the parame- 
ter goodness of fit in 3+2 schemes with respect to 3+1 
schemes is mainly a statistical effect due to an increase of 
the number of parameters, in agreement with the results 
of Ref. [l^. Therefore, in this paper we consider only 
3+1 neutrino mixing. 

We analyzed the MiniBooNE antineutrino data ex- 
tracted from the figures presented in Refs. [3l|, [s^. We 
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FIG. 2. Allowed regions in the sin^ 2'!?efi-Am4i plane ob- 
tained from the fit of MiniBooNE antineutrino data with en- 
ergy E > 475 McV 31, 32]. The best-fit point at sin^ = 
0.0045 and Am|i — 4.79 eV^ is indicated by a cross. 
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FIG. 3. Fit of MiniBooNE antineutrino data [3l|,|33| (points 
with error bars). The red solid line corresponds to the best 



fit (sin^2i9e 



0.005 and Am|i 



5eV' 



The blue 



dashed and green dotted lines correspond, respectively, to: A: 



sin^ 2i9e 



0.005 and Amii = 0.8 eV^ B: sin^ 2i9e 



0.01 



and Amii = 0.5 eV 



followed the method described in the MiniBooNE web 
page of the data release relative to Ref. [1^ , modified by 
rescaling the predicted signal by the ratio of POT. Fur- 
thermore, in order to reproduce the allowed regions in 
the sin 2i?-Arn^ plane presented in Refs. '32f|, we 
increased the background uncertainty by a factor 1.1. 
In the fit we considered also the data obtained with 
5.66 X 10^" POT [iBi, which are important because of the 
correlated uncertainties of Pe and i^^ data. 

The results of our fits are shown in Figs. [T] and [21 re- 
spectively, for all MiniBooNE antineutrino energy bins 
and for the energy bins above 475 MeV. One can see 
that the allowed regions are similar to the corresponding 
ones presented in Refs. [3ll.l3^. 

As one can see from Fig. [3l the low-energy anomaly 
is not well fitted for the best-fit value of the oscillation 
parameters, sin^ 2?9e^ = 0.005 and Am|]^ = 4.68 eV^, but 
can be fitted in the case of a lower value of Am|]^. In 
the example A (sin^ 2i?e,x = 0.005 and Amj^ = 0.8 eV^ 
which is within the 68% C.L. allowed region in Fig. [1]) 
the value of sin^ 2-(?e/j is the same as in the best-fit point 
but the lower value of Aml^ increases (Pp^^p^) in the 
low-energy bins. In the example B (sin^ 2z?ep = 0.01 
and Am|j^ — 0.5 eV^, which is also within the 68% C.L. 
allowed region in Fig. [IJ the larger value of sin^ 2?9eM 
and the smaller value of Amli allow us to fit the low- 
energy anomaly even better. Although by eye the lines 
corresponding to the cases A and B may appear to fit all 
the MiniBooNE antineutrino data better than the best- 
fit line, the best-fit point has a lower value of because 



of the correlations of the uncertainties of the bins, which 
are given by the covariance matrix of the MiniBooNE 
data release relative to Ref. [1^ . 



III. MINOS 

The MINOS collaboration presented recently [s^ the 
updated results of a search for j/^ disappearance into ster- 
ile neutrinos obtained by comparing the samples of neu- 
tral current (NC) events measured at the near detector 
(ND) and far detector (FD). In the MINOS experiment 
the neutrino beam is produced through the decay of pions 
generated by 120 GeV protons hitting a graphite target. 
The pions fly in a 675 m long decay pipe. The near and 
far detectors are located, respectively, at the distances 
^ND = 1.04km and LpD = 735km from the target. The 
analysis presented in Ref. |33|] limits |C/;i4p below 0.019 at 
90% C.L. assuming that there are no oscillations before 
the near detector and that the oscillations are completely 
averaged in the far detector. Since the neutrino energy 
range goes from about 1 GeV to about 20 GeV, the first 
condition is satisfied for Airily < 1 eV^ and the second 
condition is satisfied for Ara\-^ > 0.2 eV^. Hence, the 
range of Amli for which the bound on |C/;^4p presented 
in Ref. [s^ is limited. 

Since we consider higher values of Am|]^, we analyzed 
the data presented in Ref. (33j taking into account pos- 
sible oscillations before the near detector. From our ex- 
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FIG. 4. Fit of MINOS data extracted from Ref. [H. The 
red solid line corresponds to the best fit (sin^ 2i9en = 0.005 
and Amli = 4.68 eV^). 



FIG. 5. Raster-scan upper bound on sin'^ 2i?^^ as a function 
of Amli plane obtained from the fit of MINOS neutral current 
data 



traction of the data in Ref. [33 
shown in Fig. U for the ratio 



we obtained the values 



R = 



FD 



(1 - P. 



;nd 



(7) 



In the fit we considered a fully correlated 2.9% system- 
atic uncertainty [33|. We also took into account as fully 
correlated the uncertainties given in Ref. on the pos- 
sible effect of appearance in the far detector due to 
|C/e3p < 0.040, which is the 2010 MINOS 90% limit 
This value is compatible with the recent results of the 
T2K ^39,] and MINO S ,40.1 e xperiments which are in fa- 
vor of Ve appearance |4ll l42j . 

We calculated the oscillation probability at the near 
detector with the approximate method derived in 
Ref. [i^, which takes into account the partial decoher- 
ence of the neutrino state at the production due to the 
fact that the decay length of the parent pion and the 
length of the decay pipe are comparable with the distance 
from the target to the near detector. For completeness, 
we took into account also possible oscillations in the far 
detector due to Am^j^, but we neglected for simplicity 
possible Vg Vs transitions due to Am§2 (see Eq. (11) 
of Ref. [4J|)- This is equivalent to assuming a negligible 
value for \Usz\- 

We averaged the oscillation probabilities over the neu- 
trino flux and the neutral current cross section, taking 



into account the energy resolution of the detector: 

/d^(T 



(8) 



where / is the same integral without P{E^). We took the 
neutrino flux <^^[E^) as a function of the neutrino energy 
El, from Fig. 2.10 of Ref. [11]. The differential cross 
section has been approximated by the deep- inelastic cross 
section on a isoscalar target using the NNLO MSTW 
2008 set of parton distribution functions 46]. For the 
energy resolution function ijj-[cs{E\i, £'roco) which connects 
the hadronic energy E\^ to the reconstructed energy -Ercco 
we used a Gaussian distribution with standard deviation 
56%//^ dl. 

Figure [5] shows the upper bound on sin^ as a func- 
tion of Arn^i plane obtained with a raster scan, which can 
be compared with that obtained by the MINOS collabo- 
ration in Ref. 



One can see that for Am^j^ 



< 1 eV" 



we have the limit sin^ 2i^^^ < 0.09 at 90% C.L., with wig- 
gles due to oscillations in the far detector. Using Eq. ([B]), 
this limit corresponds to |?7;^4p ^ 0.023, which is about 
the same as that obtained by the MINOS collaboration 



in Ref. 
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For 



Am|i 



> leV^ the upper bound on 
sin"^ 2'3nn rapidly disappears, i n ag reement with the dis- 



cussion above and that in Ref. [43j. 

Figure [5] shows the exclusion curves in the sin^ 2i? 
Amli plane obtained with a two-parameters least- 
squares analysis. Also in this figure one can see that 
the oscillations in the near detector weaken the bound 



for 



Arn|;^ 



> leV . The best-fit point at 
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FIG. 6. Exclusion curves in the sin^ 2'dnfj,-Ain1i plan e ob- 



tained from the fit of MINOS neutral current data 



The 



best-fit point at sin 2i9fi;j = 0.48 and Am4i — 6.76 eV is 
indicated by a cross. 



sin^ 2^, 



0.48 and At 



'41 



= 6.76 eV^ allows to fit the 



small positive values of R measured in the low-energy 
bins, which have the smaller uncertainties. This effect is 
in agreement with the discussion in Ref. p3l |. where it 
is explained that for a value of lS.m\^ such that the low- 
energy bins correspond to the first oscillation minimum 
at the near detector the denominator in Eq. ([7]) is smaller 
than the numerator, leading to i? > 1. 



IV. GLOBAL 3+1 FITS 

In this section we present the results of the update 
of our analysis in Ref. [s^ which takes into account the 
new MiniBooNE antineutrino data discussed in Sectionllll 
and the bound on |?7;^4p obtained from MINOS neutral- 
current data in Section IIIII As in Ref. [s^, we con- 
sider also the short-baseline 't'^ — ^-Ve data of the LSND 
[H, KARMEN i23|], NOMAD [H and MiniBooNE neu- 
trino psl j experiments, the short-baseline disappear- 
ance data of the Bugey-3 HI, Bugey-4 [50|, R0VN091 
[5l| . Gosgen [s^], ILL [H^ and Krasnoyarsk |5^] reactor 
antineutrino experiments, taking into account the new 
calculation of the reactor flux [551,^561] which indicates a 
small Ue disappearance (the reactor antineutrino anomaly 



yses of all these data without and with the data of Gal- 
lium radioactive source experiments (GALLEX [60l - [6^ 
and SAGE (s^-E^) which indicate a Ve_ disappearance 
(the Gallium neutrino anomaly [isl. [stI, [67i-'7?| ) . We an- 
alyze the Gallium data according to Ref. '48] . 

For the new MiniBooNE antineutrino data discussed 
in Section In] we consider two cases: 

LOW: All MiniBooNE neutrino and antineutrino data, 
including the low-energy bins from 200 MeV to 
475 MeV. 

HIG: Only MiniBooNE neutrino and antineutrino data 
with energy E > 475 MeV. 

TableUshow the results of our global analyses for these 
two cases without and with Gallium data (GAL). The 
corresponding allowed regions in the siv? 2-df.fi-Am\i, 



1.^1 and sin 



planes are shown 



sin 2i^ee-Amj 
in Fig. [71 

From Tab. U one can see that in all cases the fit of 
the data without oscillations is rather bad and the os- 
cillations in the framework of 34-1 mixing have a good 
goodness of fit. The goodness of fit is larger in the HIG 
case than in the LOW case, because the best fit of Mini- 
BooNE data does not fit well the low-energy bins, as we 
have seen for the antineutrino data in Section|lIl Gallium 
data slightly worsen the goodness of fit because the best 
fit of the analysis of Gallium data [1^ requires a value 
of I C/e4 1 ^ which is larger than the bound given by reactor 
antineutrino data (see Ref. [SO]*). 

The best-fit value of the 3-1-1 oscillation parameters 
is rather stable under variations of the considered data 
sets It points to a rather large value of Amli, which 
in a hierarchical scheme corresponds to « 2.4 eV. 
This value is in tension with the limits given by standard 
ACDM analyses of cosmological data _14- 19i]. If it will be 
confirmed by future data, it may indicate the existence of 
non-standard effects in the evolution of the Universe [l| . 
However, as one can see from Fig. [7] there are allowed 
regions at Amj^i ss 1 eV^ which are more compatible with 
standard ACDM cosmology. For example, in the HIG fit 

'41 ■ ~ ^ 

~ 0.0012, sin2 2i9e ~ ""^ "-"^ ' 
ill ~ 1.2 eV^ sin^ 2de^ - " ""^ ' "-"^ ' 
sin2 2i?„„ w 0.051; C) Am|, w 0.91 eV^ sin^ 2i9, 

,2 < 



there are three regions within la at: A) Amli « 1.6 eV^ 

0.12, sin2 2i9^^ « 0.037; 



sin^ 2-de 
B) Ami 



0.0014, sin^ 2de 
V 0.91 eV^ 
0.078. 



0.11, 



[13), the KamLAND ^ bound on |C/e4p (see Ref. (Hj), small values of sin^ 2^, 



the short-baseline disappearance data of the CDHSW 
experiment js^ and the constraints on |J7;i4p obtained 
in Ref. [s^l from the analysis of the data of atmospheric 
neutrino oscillation experiments. We present global anal- 



0.0020, sin^ 2'&ee « 0.10, sin^ 2i?^^ 

Figure [7] shows also the 3cr contours of the regions al- 
lowed by appearance (APP) and disappearance (DIS) 
data. One can see that they are compatible with the 
global allowed regions. It is interesting to note that 
the values of sin^ 2??ee and sin^ 2??^^ allowed by the 
analysis of appearance data can be much smaller than 
those allowed by the global analysis. The reason is that 

and sin^ 2i?^^ can be obtained 
not only with small values of |C/e4p and |f/;i4p, respec- 
tively, but also with |t/e4p and |C/^i4p close to unity 
(see Eq. dH])). Since one can fit the appearance data 



with \Ue4 



1 and sin 2i}p 
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LOW LOW+GAL HIG HIG+GAL 

No Osc. 174.8 186.8 157.8 169.8 

NDF 133 137 127 131 

GoF 0.0088 0.003 0.033 0.013 



3+1 


2 

Xmin 


134.9 


142.2 


120.7 


128.0 




NDF 


130 


134 


124 


128 




GoF 


0.37 


0.30 


0.57 


0.48 




AmJi [eV J 


5.6 


5.6 


5.6 


5.6 




0.032 


0.037 


0.033 


0.038 




sin^ 2i9e,i 


0.014 


0.012 


0.013 


0.011 




0.0018 


0.0018 


0.0017 


0.0017 




sin^ 2l9ee 


0.12 


0.14 


0.13 


0.14 




sin^ 2i?^^ 


0.054 


0.049 


0.05 


0.045 


PG 


AXmin 


15.8 


15.8 


9.3 


9.2 




NDF 


2 


2 


2 


2 




GoF 


4 X 10"" 


4 X 10"'' 


0.01 


0.01 



TABLE L Values of x^, number of degrees of freedom (NDF), goodness-of-fit (GoF) and best-fit values of the 3+1 oscillation 
parameters obtained from global fits of the data of short-baseline neutrino oscillation experiments. The first three lines 
correspond to the case of no oscillations (No Osc). The following nine lines correspond to the case of 3+1 mixing. The 
last three lines give the parameter goodness-of-fit (PG) f47| obtained by comparing the global best-fit with the best fits of 
the appearance and disappearance data. LOW and HIG refer, respectively, to MiniBooNE data with and without the three 
low-energy bins from 200 MeV to 475 MeV. GAL refers to Gallium radioactive source experiment data analyzed according to 
Ref. gg. 



and siT?2depL — 4|C/e4p, small values of sin^ 2i?ee and 
sin^ 2^^^ are allowed. On the other hand, in the global 
analysis values of |C/e4p and |t^^4p close to unity are 
forbidden, respectively, by the observation of solar and 
atmospheric neutrino oscillations. In our analysis we use 
the data of the very- long-baseline KamLAND reactor an- 
tineutrino experiment 15^ which measured a disappear- 
ance of Ve due to Ar7i2]^ which is compatible with solar 
neutrino oscillations. The KamLAND measurements re- 
quire a relatively large value of \Uei\^ + |f/e2p, which by 
unitarity constrain |C/e4p to be much smaller than unity 
(see Ref. [11]). In a similar way, the observation of atmo- 
spheric oscillations due to Amg^ requires a relatively 
large values of |C/^ip + |/7^2p + W^^?, which by uni- 
tarity constrain |C/;i4p to be much smaller than unity. 
In our analysis we use the bound on |J7^4p obtained in 
the analysis of atmospheric neutrino data presented in 
Ref. '37^. 

In Fig. [7] one can see that for lS.rri\i > 1 eV^ large val- 
ues of sin^ 2i9^^ are not allowed by the analysis of appear- 
ance data. The reason is that, as explained in Section HIl 

we fit the MiniBooNE 'i^^ -> data together with the 
MiniBooNE — data which have correlated uncer- 
tainties. The MiniBooNE — > data constrain the 

disappearance ofV^'s for Amli > 1 eV^ [ill, limiting 
the allowed value of sin^ 2i?^^. 

The last line of Tab. H] gives the parameter goodness- 
of-fit [l^l in the four analyses, which has been obtained 
by comparing the global best-fit with the sum of the 
best fits of the appearance and disappearance data. We 
think that this parameter goodness-of-fit is more reliable 



than the parameter goodness-of-fit obtained by compar- 
ing data in favor and against short-baseline neutrino os- 
cillations, which has beenj;onsidered in several previous 
analyses, including ours j30| . The reason is that the dis- 
tinction between appearance and disappearance data is 
made a priori, without considering the data. In this case 

^Xmin ~ (Xmin)app+dis — (Xmin)app ~ (Xmin)dis IS & ran- 
dom variable with a distribution with 2 degrees of 
freedom, as shown in Ref. [l^] . On the other hand, com- 
paring data in favor and against short-baseline neutrino 
oscillations one chooses the two sets of data in order to 
obtain always the worse Ax^in allowed by the data. In 
this case Ax^i^ does not have a distribution and the 
parameter goodness-of-fit is underestimated by assuming 
a x^ distribution with 2 degrees of freedom. Moreover, 
after the discovery of the reactor antineutrino anomaly 
(S?! it is not clear if the reactor antineutrino data should 
be put in the group of experiments in which there is no 
evidence of short-baseline neutrino oscillations, as tradi- 
tionally done, or in the group of experiments in favor of 
short-baseline neutrino oscillations. Therefore, we advo- 
cate the robust parameter goodness-of-fit [IJ obtained 
from appearance and disappearance data. 

The values of such parameter goodness-of-fit listed in 
the last line of Tab. U show that there is a tension be- 
tween appearance and disappearance data. This tension 
is severe in the two LOW fits. The reason can be un- 
derstood by looking at the two corresponding panels in 
Fig. [71 One can see that for appearance data there is no 
3ct allowed region around the best-fit point, because the 
low-energy MiniBooNE data are not fitted well by such 
high values of Aml^ and small mixing. Hence the best 
fit point lies out of the region of overlap of the 3a allowed 
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FIG. 7. Allowed regions in the sin^ 2i9e^-Am4i, sin^ 2i9ee-Am|i and sin^ 2i9^^-Am4i planes in the four cases listed in Table|Tl 
The best-fit points are indicated by crosses. The thick solid blue lines with the label APP show the 3cr allowed regions obtained 

fi-om the analysis of*!//l -^''fl appearance data. The thick solid red lines with the label DIS show the 3a allowed regions obtained 
from the analysis of disappearance data. 
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regions of appearance and disappearance data. This is a 
symptom of a severe tension. 

The tension between appearance and disappearance 
data is reduced in the HIG fits, for which there is a 3cr 
allowed region around the best-fit point, as one can see 
from the corresponding panels in Fig. [T] The resulting 
parameter goodness-of-fit, about 1%, is not large, but 
also not small enough to reject the fit with reasonable 
confidence. 

Table |T] and Fig. [7] show that the Gallium data do not 
have a large impact on the results of the fit. The reason 
is that the data points are only four, much less than the 
reactor data points which give information on the same 
probability of and Ve disappearance. The main effect 
of Gallium data is to favor values of Amli larger than 
about leV^. Indeed, in the HIG+GAL panel in Fig. [7] 
there are no regions allowed at la below about 1.3 eV^. 

In comparison with the allowed regions of the 3+1 os- 
cillation parameters obtained in Ref. ^30J , the inclusion in 
the analysis of the MINOS data discussed in Section Hill 
has the effect of disfavoring the regions with Amli < 
1 eV^ and moving the best-fit value of Am|j^ from the 
value of 0.9eV^ obtained in Ref. [sO*] to Am^i w 5.6 eV^ 



V. CONCLUSIONS 

The results of our analysis of short-baseline neutrino 
oscillation data show that the data can be fitted in the 
framework of 3+1 neutrino mixing, which requires the 
existence of a sterile neutrino with mass at the eV scale. 

The fit has a tension due to the lack of observation 
of enough short-baseline disappearance of Dg and to 
fully explain the — >■ Pe signal observed in the LSND 
and MiniBooNE experiments. However, we found that 
the appearance and disappearance data are marginally 
compatible if we neglect the data on the MiniBooNE 
low-energy anomaly, which may have an explanation dif- 
ferent from'i/^ — ^-'j^e oscillations. In any case, we think 
that the neutrino oscillation explanation of the data can- 
not be dismissed with a light heart, because besides the 
LSND and MiniBooNE indications in favor of a short- 
baseline Du, signal we have the reactor antineutrino 
anomaly [57[ and the Gallium neutrino anomaly |48jJ in 
favor, respectively, of short-baseline v^, and v^^ disappear- 
ance which could be due to the same squared-mass dif- 



ference. 

Since the recent MiniBooNE antineutrino data [3l|, H^l 
do not show a large difference from the neutrino data 
[25} . there is no serious motivation to consider the more 
complicated 3+2 neutrino mixing, which would allow for 
a possible CP-violating difference between neutrino and 
antineutrino transitions Wk [29l . Isol . l33 - [37j . Moreover, 
as we have shown in Ref. 3C|. 3+2 mixing cannot re- 
solve the tension between appearance and disappearance 
data. Finally, the hierarchical 3+1 scheme (see Eq. ([T])) 
is favored over a hierarchical 3+2 scheme by standard 
ACDM analyses of cosmological data p^ - [T9| , which dis- 
favor sums of neutrino masses much larger than 1 eV, 
and by Big-Bang Nucleosynthesis data, which allow the 
existence of one sterile neutrino [2^ Hll, but not more 
[19I [22I (keeping however in mind the caveat that these 
bounds refer to the number of sterile neutrinos which 
have been fully thermalized in the early Universe). 

Hence, in this paper we considered only 3+1 neutrino 
mixing, which is attractive for the natural correspon- 
dence of the existence of one new entity (a sterile neu- 
trino) with the observation of a new effect (short-baseline 
oscillations). 

The results of our fit excluding the MiniBooNE low- 
energy anomaly lead to a best fit at Am^i ~ 5.6 eV^, 
which is larger than that obtained in Ref. [s^ (about 
0.9 eV^) because of the new MINOS constraints dis- 
cussed in Section IIIII The new best fit value of Am|j^ 
is rather large in c omp arison with the standard cosmo- 
logical bounds [lj-[l9| and may indicate the existence 
of non-standard effects in the evolution of the Universe 
[igj . However, there are three regions within la at 
Amli « 1.6, 1.2, 0.91 eV^ which may be compatible 
with the standard cosmological bounds. 

We have also shown that the data on the Gallium neu- 
trino anomaly favor values of Am\^ larger than about 
1 eV^ , but their impact is small because the results of 
the analysis are dominated by the more abundant re- 
actor antineutrino data, which give information on the 
same probability of v,, and disappearance. 
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